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ABSTRACT 

We study the evolution of a low mass x-ray binary coupling a binary stellar evolution 
code with a general relativistic code that describes the behavior of the neutron star. 
We assume the neutron star to be low-magnetized (B ~ 10 8 G). In the systems 
investigated in this paper, our computations show that during the binary evolution 
the companion transfers as much as 1 Mq to the neutron star, with an accretion rate of 
~ 10 -9 M Q /yr. This is sufficient to keep the inner rim of the accretion disc in contact 
with the neutron star surface, thus preventing the onset of a propeller phase capable 
I of ejecting a significant fraction of the matter transferred by the companion. In this 

. scenario we find that, for neutron stars governed by equations of state from soft up to 

moderately stiff, an accretion induced collapse to a black hole is almost unavoidable. 
The collapse to a black hole can occur either during the accretion phase or after the 
end of the mass transfer when the neutron star is left in a supramassive sequence. In 
this last case the collapse is driven by energy losses of the fast spinning magneto-dipolc 
rotator (pulsar). For extremely supramassive neutron stars these energy losses cause 
a spin up. As a consequence the pulsar will have a much shorter lifetime than that of 
a canonical, spinning down radio pulsar. This complex behavior strongly depends on 
the equation of state for ultra-dense matter and therefore could be used to constrain 
the internal structure of the neutron star. In the hypothesis that the r-modes of the 
r> I neutron star are excited during the accretion process, the gravitational waves emisson 

^ ■ limits the maximum spin attainable by a NS to roughly 2 ms. In this case, if the mass 

transfer is conservative, the collapse to a black hole during the accretion phase is even 
more common since the maximum mass achievable before the collapse to a black hole 
during accretion is smaller due to the limited spin frequency. 

Key words: Stars: neutron - X-rays: binaries - binaries: close - pulsars:general - 
relativity 
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1 INTRODUCTION 

The widely accepted scenario for the formation of millisec- 
ond radio pulsars (MSP) is the recycling of an old neu- 
tron star (hereafter NS) by a spin-up process. The spin-up 
is due to the accretion of matter and angular momentum 
from a Keplerian disc that is fueled via Roche lobe over- 
flow of a binary late-type companion (see Bhattacharya & 
van den Heuvel 1991 for a review). If enough mass and an- 
gular momentum are transferred, the NS spin attains an 
equilibrium value that is roughly equal to the keplerian an- 
gular frequency at the inner rim of the accretion disc (Ghosh 
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& Lamb, 1979). Since the NS has a weak surface magnetic 
field (~ 10 8 G), the magnetospheric radius (at which the disc 
pressure is balanced by the magnetic pressure) truncates the 
accretion disc close by or at the NS surface, and the equilib- 
rium period is expected to be, in most cases, below one mil- 
lisecond. Typically ~ 0.35 Mq are sufficient to reach a spin 
period of 1 ms (e.g. Burderi et al. 1999). Most donor stars 
in systems hosting recycled MSPs have certainly lost, dur- 
ing their interacting binary evolution, a mass greater than 
0.35 A^q since they now appear as white dwarfs of mass 
0.15 - 0.30 M Q (e.g. Taam, King, & Ritter 2000), whose 
progenitors are likely to have been stars of 1.0 — 2.0 Mq 
(Webbink, Rappaport, & Savonije 1983; Burderi, King, & 
Wynn 1996; Tauris & Savonije 1999). Therefore if the mass 
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transfer is conservative the amount of matter accreted is 
well sufficient to spin the star up to periods below 1 ms 
(Cook, Shapiro, & Teukolsky 1994a), or even to produce an 
accretion induced collapse into a black hole. Once the accre- 
tion and spin-up process ends, the magnetospheric radius 
expands beyond the light cylinder radius (where an object 
corotating with the NS attains the speed of light). This ini- 
tiates a phase in which the rotational energy of the NS is 
emitted via electromagnetic radiation and the star can be 
observable as a rapidly rotating radio pulsar. According to 
this model, Low Mass X-ray Binaries (hereafter LMXBs) 
are the progenitors of MSPs. Indeed the discovery of coher- 
ent X-ray pulsations in four transient LMXBs, namely SAX 
J1808. 4-3658 with a spin period P — 2.5 ms and an orbital 
period of P orb = 2 h (Wijnands & van der Klis 1998), XTE 
J1751-305 (P = 2.3 ms, P orb = 42 min, Markwardt et al. 
2002), XTE J0929-314 (P = 5.4 ms, P orb = 43 min, Gal- 
loway et al. 2002), and XTE J 1807-294 with a spin period 
P = 5.2 ms and an orbital period of P or b = 40 min (Mark- 
wardt et al. 2003) confirmed that NSs in LMXBs can be 
accelerated to millisecond periods. 

Several numerical methods have been developed to solve 
the Einstein equations for a rotating NS (see Stergioulas 
1998 for a review). Stability criteria show that a rapidly 
rotating NS can support a maximum mass (against gravi- 
tational collapse) much higher than the non-rotating mass 
limit, since the centrifugal force attenuates the effects of 
the gravitational pull (e.g. Friedman, Ipser, & Sorkin 1988). 
Conversely, if a rotating NS has a mass that exceeds the non- 
rotating limit (i.e. a supramassive NS), it will be subject to 
gravitational collapse if it loses enough rotational energy. 
Numerical simulations of rotating NSs show that, contrarily 
to the standard behavior, supramassive NSs spin up just be- 
fore collapse, even if they lose energy. This effect is known 
to be stronger for higher mass objects (Cook, Shapiro, & 
Teukolsky 1992). 

The value of the maximum rotating and non-rotating 
mass depends on the equation of state (EOS) governing the 
NS matter. On the other hand, the minimum allowed pe- 
riod for a given mass occurs when gravity is balanced by 
centrifugal forces at the NS equator (mass shedding limit). 
Thus the spin period can be used to constrain the mass- 
radius relation for the NS, i.e. its EOS. In the context of the 
standard (gravitationally-bound) NS models (e.g. Glenden- 
ning 2000), several EOSs have been proposed. We usually 
distinguish different EOSs depending on their stiffness (i.e. 
the value of dp/de, were p is the fluid pressure and e is the 
energy density) . If the EOS is stiff, the matter is less com- 
pressible at high densities, resulting in larger NS radii as 
compared to a soft EOS, and hence in longer minimum ro- 
tational periods. Except for few, very stiff cases, most EOS 
predict minimum rotational periods well below 1 ms. How- 
ever, no sub-millisecond pulsars have been detected up to 
date: the shortest observed spin period is ~ 1.5 ms (Backer 
et al. 1982), uncomfortably higher than the theoretical pre- 
dictions. 

In an attempt to find an explanation for the apparent 
clustering of the spin periods of millisecond pulsars around 
2 ms, Bildsten (1998) and Andersson (1998) indipendently 
suggested that LMXBs emit gravitational waves once they 
reach a critical spin frequency. Burderi and D'Amico (1997) 
showed that for nonaxisymmetric m-modes, assuming a re- 



alistic range of temperatures, the values of the critical spin 
frequency are remarkably close to the limiting spin frequency 
determined by the centrifugal limit at the border of the NS. 
On the other hand, Andersson, Kokkotas and Stergioulas 
(1999) demonstrated that at a certain spin frequency (much 
lower than the maximum attainable spin period) an insta- 
bility to the Rossby waves (r-modes) of the star arises, thus 
causing emission of gravitational waves. Levin (1999) sug- 
gested that the gravitational waves emission causes the on- 
set of a spin-up spin-down cycle, and not of a steady state 
spin equilibrium: in this scenario the NS undergoes a very 
rapid spindown (lasting ~ 1 yr) due to the rapid heating 
during the r-mode excitation, and then starts another cycle 
of accretion driven spin up. Brown and Ushomirsky (2000) 
showed that if the NS has a superfluid core the steady state 
scenario is not viable because the predicted quiescence lumi- 
nosity in this case is much higher than the observed one in 
the X-Ray transient Aql X-l. Levin and Ushomirsky (2001) 
showed that, when keeping into account the presence of the 
solid crust, the critical spin frequency for the onset of the 
r-mode varies between ~ 600 Hz and ~ 200 Hz, depending 
on the core temperature. 

In this paper we consider the full evolution of a LMXB, 
trying to determine how the results of our modeling of the 
recycling scenario compare to the observations and which 
effects peculiar to General Relativity are indeed observable. 
We will also show the differences in the evolution of the sys- 
tem when the r-mode instability is excited during accretion 
and when it is not excited. 



2 EVOLUTION EQUATIONS FOR THE 
COMPACT OBJECT 

A rotating NS is unambiguously defined by the boundary 
conditions for the integration of the Einstein equations or, 
equivalently, by suitable pairs of resulting integrated quanti- 
ties, such as the baryonic mass and the angular momentum 
or the baryonic mass and the total mass-energy of the star. 
Therefore the evolution of the NS is determined by the evo- 
lution of such pairs. We consider the evolution of the NS 
both during the accretion phase and after the accretion has 
finished: in the former case it is convenient to solve the evo- 
lution equations for baryonic mass and angular momentum, 
since we know the general formula for the torque exterted by 
the accreting matter on the NS, while in the latter we solve 
the evolution equations for baryonic mass and mass-energy, 
since we can evaluate the luminosity of a magnetodipole ro- 
tator using Larmor's formula. 

2.1 Evolution equations for the accretion phase 

According to accretion theories, matter transferred from the 
companion star to the NS via Roche- Lobe overflow forms an 
accretion disc around the compact object. The disc is trun- 
cated because of one of the following reasons: i) the inter- 
action with the magnetic field of the NS, which truncates 
the disc at the magnetospheric radius r m ; ii) the presence 
of NS surface itself at Pins; iii) the lack of closed Keplerian 
orbits for radii smaller than the marginally stable orbit ra- 
dius, Pmso (at few - depending on the mass and spin of the 
compact object - Schwarzschild radii from the NS centre). 
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The magneto-spheric radius is defined as the radius at 
which the pressure of the disc equals the pressure of the 
magnetic field of the NS. The magnetospheric radius can 
be written as a fraction <j> (see Burderi et al. 1998) of the 
Alfven radius _Ra (the radius at which the energy density of 
the, assumed dipolar, NS magnetic field equals the kinetic 
energy density of the spherically accreting matter) as 
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where a < 1 is the Shakura-Sunyaev parameterization of the 
accretion-disc viscosity (for which we will assume a typical 
value of 0.1), ma is the NS gravitational mass in Mq, rhs 
is the baryonic accretion rate on to the NS in 10 -8 Mq /yr, 
and Ra is 



Ra = 1.24 x lQ a ii% 7 m~ in m~ 2/7 cm 



(2) 



where /i26 is the magnetic dipole moment of the NS in units 
of 10 26 G cm 3 , defined from B s = fi/R 3 , where R is the 
equatorial radius and B s the surface magnetic field of the 
NS at its equator. 

The NS radius, that is always in the order of 10 6 cm, 
depends both on the mass of the NS and on its angular 
momentum. In general smaller radii correspond to larger 
masses while larger radii correspond to larger angular mo- 
menta. Thus a rapidly rotating NS can have a much larger 
radius than a non-rotating one (the equatorial radius can 
expand up to 40%, see Cook, Shapiro, & Teukolsky 1994b). 

The marginally stable orbit is the smallest stable orbit 
possible for a test particle around an axisymmetric, rotating 
body of gravitational mass M G and angular momentum J. 
Its radius, following Bardeen (1970), is: 



Rmso 
Z 2 



= R g {3 + Z 2 - yj(3 - Zi)(3 + Zi + 2Z 2 )); 



Z x = 



R- 




where c is the speed of light and G is the gravitational con- 
stant. Therefore the inner radius of the disc, Rd, will be: 



Rd = max {r m , .Rns, -Rmso} • 



(3) 



The position of the inner rim of the disc is crucial for the 
accretion of matter on to the NS. In particular it is very im- 
portant in the case Rd = r m : since at r m the matter is forced 
by the magnetic field to corotate with the NS, accretion on to 
a spinning magnetized NS is centrifugally inhibited once r m 
lies outside the corotation radius r co , the radius at which the 
Keplerian angular frequency of the orbiting matter is equal 

fi 1/3 2/3 

to the NS angular frequency: r co = 1.50 x 10 m G P_ 3 cm, 
where P_3 is the spin period in milliseconds. Conversely if 
the magnetospheric radius is smaller than the corotation 
radius accretion of matter can proceed undisturbed. How- 
ever, as r m scales as a negative power of mj, a decrease 



in the mass transfer rate, which can occur for instance at 
the end of the accretion phase, will results in an expansion 
of the magnetosphere. In this paper we considered NSs with 
/^26 — 1, which is typical of observed MSPs (see e.g. Lorimer 
1994), corresponding to surface magnetic fields of ~ 10 8 G. 
Moreover, in the cases considered here, the accretion rate is 
always above ~ 10 _10 Mo/yr during the mass transfer phase 
(see figure 0. In this case, eq. Q implies a magnetospheric 
radius of ~ 10 6 cm, that is comparable with the NS radius. 
Although we considered in our simulations (see below) the 
possibility of a propeller phase, in these conditions the de- 
tails of the interaction of the accretion disc matter with the 
NS magnetosphere during the accretion phase will not affect 
much our results. In fact the magnetospheric radius lies out- 
side of the corotation radius only for a small time ( <0.1% 
of the evolution time) just before of the end of accretion (see 
figure [5] for details). In other words, the spin evolution of a 
low magnetized NS is not very different from the spin evolu- 
tion of a non-magnetized NS, in which no propeller effect is 
possible since the disc is always truncated near the surface 
of the NS. We can therefore write the mass evolution equa- 
tion for r m < r co , which is simply given by the condition 
that the baryonic mass per unit time accreted on to the NS 
is equal to the mass lost by the companion per unit time 
due to Roche-Lobe overflow, Mrl, i.e. 



M B = M E 



(4) 



Lamb, Pethick, & Pines (1973) wrote a general equa- 
tion describing the flow of angular momentum into the stel- 
lar magnetosphere including the material stress at the inner 
edge of the disc as well as magnetic and viscous stresses in a 
disc partially threaded by the NS magnetic field. In partic- 
ular they assumed that the threading occurs in a transition 
region near the inner rim of the disc where the magnetic 
field of the NS couples to the accretion disc. However they 
noticed that, for slow rotators, magnetic and viscous stresses 
can be neglected with respect to the contribution of the ma- 
terial stress. Ghosh, Lamb, & Pethick (1977) showed that, 
in the case of a disc rotating in the same sense of the star, if 
the transition region where the magnetic field couples to the 
disc is small, these additional contributions are still negligi- 
ble even for rapidly rotating NS. In the present discussion 
we assume therefore that the torque exerted on the NS is 
only due to the contribution of the material stress at the 
inner edge of the disc (Pringle & Rees 1972). We will study 
the evolution of high-magnetized NSs with heavily threaded 
accretion discs in a future paper. Note, however, that Wang 
(1996) showed that threading effects can modify the max- 
imum achievable period (i.e. the equilibrium period) of the 
NS only by a few percents with respect to the unthreaded 
case. 

The angular momentum per unit baryonic mass of a 
particle orbiting around a rotating, axisymmetric object at 
a distance R is (see Bardeen 1970) 
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Thus we can write the equation for angular momentum evo- 
lution of the NS simply as 



J = j(R D )-M B , 
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Figure 1. Baryonic mass accretion rates (in MQ/yr) as a function of time (in yr) for systems consisting of a NS with initial gravitational 
mass of 1.4Mq and: a) a population II (low metallicity) star of 0.85Mq transferring about 0.65Mq; b) a 1.15Mq population I star 

ring about O.99M . 
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Figure 2. Time evolution of the relevant radii for the accretion 
process on to the neutron star in the binary system of figure lb. 
The EOS governing the ultradense matter is EOS FPS. The solid 
line is the corotation radius, the thick dashed line is NS radius, 
the dot dashed line is the magnetosphcric radius, the dotted line 
is the radius of the marginally stable orbit. Note that the mag- 
netosphcric radius exceeds the corotation radius only at the end 
of accretion, when the mass transfer rate decreases, and that for 
the first part of the accretion process the last stable orbit is out- 
side the NS. The magnetosphcric radius data are smoothed in 
order to avoid the disturbing visual effects due to the numerical 
instabilities of Mg . 



The first equation does not apply when the NS is at the mass 
shedding, i.e. the regime in which the matter at the border 
of the star has the Keplerian velocity at that radius; for each 
value of baryonic mass, the mass shedding regime is individ- 
uated by the corresponding maximum angular momentum 
of the star, JmaxiMs)- At mass shedding, the matter of 
the disc should dissipate angular momentum to accrete on 
to the star. If the disc matter has an angular momentum 
that will make the star exceed the mass shedding limit, it 
will just not be able to accrete since it will not be gravi- 
tationally bound to the star, until viscous forces drive the 
excess of angular momentum to the outer zones of the disc 
(where it will be given back to the companion through tidal 
forces) allowing matter to accrete on to the star. In this 
situation, the second equation holds. Integrating it we get 
Jf — Ji = Jmax,/ - Jmax.i, but since J; = Jmax.i we have 
Jf = Jmaxj- the accreted matter will give the NS only the 
angular momentum that keeps it at mass shedding, and the 
star will continue to move along the maximum rotation line. 



2.2 Evolution equations at the end of the 
accretion phase 

The NS is thought to switch on as a radio pulsar when the 
inner edge of the disc lies outside the light cylinder radius 
(i.e. the radius at which a particle corotating with the star 
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will have velocity c), r; c = c/luns - where lons is the an- 
gular velocity of the NS: this certainly happens when the 
accretion stops and thus the disc disappears. The emission 
mechanisms for a radio pulsar are believed to be rotating 
magnetic dipole radiation and magnetospheric currents as- 
sociated with the emission of relativistic particles, both de- 
pending on the angle i between the NS magnetic moment 
and the spin axis (Goldreich & Julian 1969). These two 
emission mechanisms compensate in such a way that the to- 
tal energy emitted is nearly independent of i (Bhattacharya 
& van den Heuvel 1991). The energy loss per unit time will 
be E = — ^jfi 2 u)% s . It is then convenient to describe the NS 
evolution in terms of baryonic mass and total energy rather 
than in terms of baryonic mass and angular momentum as 
we did for the accretion phase. We can write the evolution 
equations as 

Mb = 

M G = -^fi 2 uj% s . (7) 

Depending on the different conditions, we use either equa- 
tions @ and @ or equations Q to compute the spin evo- 
lution of the NS. 

2.3 Excitation of r- modes 

Since our main goal is to compare qualitatively the evolution 
of systems in which the gravitational waves damping is not 
present and systems in which it is present, we will not go 
into the details of the theory of r-modes excitation. We will 
simply assume that, if the r-mode mechanism is present, any 
star attaining a period of 2 ms (corresponding to a frequency 
of 500 Hz) undergoes a rapid spindown to a period of 5 ms 
due to gravitational waves emission following the scenario 
proposed by Levin (1999). After this rapid spin-down, the 
NS starts again its accretion driven spin up process. 

2.4 Numerical methods of integration 

During accretion phases we coupled through equations Q 
and a detailed description of the binary evolution of the 
system obtained with the stellar evolution code with a fully 
relativistic calculation of NS physical properties given its 
EOS, its baryonic mass and its angular momentum. 

The evolution of the binary system is followed selfcon- 
sistently including the full computation of the structure of 
the secondary star, by means of the ATON1.2 stellar evolu- 
tion code (D'Antona, Mazzitelli & Ritter, 1989). The equa- 
tions of stellar structure are numerically solved by a full 
Newton-Raphson integration from the center up to the ba- 
sis of the stellar atmosphere. The numerical inputs are de- 
scribed in Mazzitelli (1989). The secondary star is consid- 
ered to be the mass losing component of the binary system, 
and its mass loss rate is computed following the formula- 
tion by Ritter (1988), as an explicit exponential function of 
the distance of the stellar radius to the Roche lobe, in units 
of the pressure scale height. The evolution of the binary 
parameters can be followed by considering several possible 
cases for the transfer of mass, and for the loss of mass and 
associated angular momentum from the system. The orbital 
evolution also includes losses of orbital angular momentum 
through magnetic braking, in the Verbunt and Zwaan (1981) 



formulation, in which the braking parameter is set to / = 1 
and through gravitational radiation. 

The relativistic computations for the compact object 
are implemented using a slightly modified version of RNS 
(Rotating Neutron Stars) public domain code by Stergioulas 
& Morsink (1999). The RNS routines provide the numeri- 
cal solution of Einstein equations for a rotating axisymmet- 
ric body, integrating the equations via Komatsu-Eriguchi- 
Hachisu method (Komatsu et al. 1989; see Stergioulas & 
Friedman 1995 for a comparison with other integration 
methods). The main problem of this approach is that we 
have evolution equations for baryonic mass and angular mo- 
mentum of the compact object, while the boundary condi- 
tions for the solution of Einstein equations are the central 
energy density and the oblateness of the star. To solve this 
problem, we used a grid of relativistic equilibrium NS mod- 
els integrating the Einstein equations for a wide range of 
initial conditions, spanning all allowed values for stable con- 
figurations. This was obtained modifying the code to build a 
complete grid of relativistic equilibrium configurations with 
the necessary numerical precision. To this we added a sta- 
bility control routine to exclude unstable configurations. 

We obtain the accretion rate Mb and a correspond- 
ing time interval (in which Mb remains unchanged) from 
the stellar evolution code. Then we integrate the differen- 
tial equations Q and JSJ using a finite differences method. 
For each integration time step the accreted baryonic mass 
is AMb = MB(t)At and the accreted angular momen- 
tum is AJ = j(R D (t))M B (t)At. For the n-th time step 
of the evolution, we search on the grid the single equilib- 
rium configuration with the updated values of Ms(t„) = 
Msitn-i) + AMb and J{t„) = J(t n -i) + AJ; in this way 
we get the corresponding values of gravitational mass, ra- 
dius, spin frequency and momentum of inertia of the NS. 
We fine-tuned the time step so that AMb, A J are always 
comparable with the distance of two neighbour points in the 
grid. 

Obviously the integration over a grid may introduce nu- 
merical uncertainties in our results: in particular, we are im- 
plicitly assuming that the spin frequency, the gravitational 
mass and the radius of the NS remain unchanged between 
two points on the grid. Therefore we introduce an error in 
the evaluation of these quantities that is equal to half the 
distance between two points on the grid, that is always well 
below 0.5% in our simulations. This also yields an uncer- 
tainty in evaluating, step by step, the value of the accretion 
radius given by eq. and the angular momentum of the 
matter at the inner rim of the disc, that is always in the 
order of a few percent. 

During the pulsar phase, we integrated equations Q in 
a similar way using the same grid of relativistic equilibrium 
NS models. To avoid the grid imprecisions we made a cubic 
spline interpolation on the values of the function uns(Mg) 
for the integration of the equation set 



3 RESULTS 

We consider three simple examples that show different pos- 
sible fates for a LMXB, depending on the characteristics 
of the companion star. In all cases we start our evolution 
with a slowly rotating (P ~ 1 s) NS with a gravitational 
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mass of 1.4 Mq - which is a typical value for a newborn 
NS (Thorsett and Chackrabarty, 1999). The EOS adopted 
to describe the ultradense matter is that proposed by Fried- 
man, Pandharipande and Smith (FPS, see e.g. Lorenz et al. 
1993 for a recent discussion) which has been widely used in 
the literature and has average values of stiffness and physical 
parameters when compared to other EOS (see Arnett and 
Bowers, 1977 for a catalog). For our simulations we consider 
the maximum mass configuration equal to the maximum ro- 
tation configuration (although this is not precisely true, see 
Stergioulas & Friedman 1995, these configurations differ one 
from each other for less than 0.1% in our case), and we have 
M max = 2.123 Mq, w max /27r = 1882 Hz. Finally, as already 
mentioned, we assume that the NS has a magnetic dipole 
moment of 10 26 G cm 3 . 

In particular we consider three typical examples of bi- 
nary evolution in which the accretion rate is low enough to 
remain conservative, both above and below the bifurcation 
period - i.e. the period below which the orbital evolution, 
during the accretion phase, proceeds towards shorter binary 
periods. The bifurcation period is P or b,bif — 18 hr for a bi- 
nary composed of a 1.4 Mq NS plus a 1.0 Mq secondary 
(Podsiadlowski et al. 2002). In the cases we considered, the 
secondary is: 

(i) A population II donor with an initial mass of 0.85A/q 
and an orbital period above the orbital bifurcation period; it 
loses 0.646Mq during the accretion phase. The initial binary 
period is 14.3 hours and it steadily increases to about 124 
h at the end of the accretion phase, which lasts ~ 1.7 x 10 9 
yr (see figure la), leaving a helium white dwarf remnant. 
Note that the orbital bifurcation period for a system com- 
posed of a 1.4 Mq NS plus a 0.85Af Q is smaller than 14.3 
hours, because the bifurcation period depends strongly on 
the companion mass. In this case, Roche-Lobe overflow is 
driven by the nuclear evolution of the companion. 

(ii) A population I donor (initially below the bifurcation 
period) with an initial mass of 1.15Mq, which loses 0.91Afo 
during the accretion phase. The orbital period evolves from 
10.6 to 3.5 hours. In this case the companion star overflows 
its Roche-Lobe due to angular momentum losses caused 
magnetic braking. The accretion phase ends when the com- 
panion becomes fully convective, according to the classical 
scenario for the evolution of cataclysmic variables in the pe- 
riod gap (Verbunt & Zwaan 1981). The accretion phase lasts 
~ 1.7 x 10 9 yr, see figure lb. 

(iii) A population I donor (initially above the bifurcation 
period) with an initial mass of 1.199M0; it loses 0.99A#q 
during the accretion phase. The period evolves from 19 to 
92 hours, leaving again a helium white dwarf remnant. The 
accretion phase lasts ~ 2.2 x 10 9 yr (see figure lc). In this 
case, Roche-Lobe overflow is driven by the nuclear evolution 
of the secondary star. 

In these evolutions the mass trasfer rate is smaller than 
the Eddington limit. Therefore it is reasonable to assume 
that the mass transfer is conservative. Moreover, the matter 
flow from the companion is not subject to large fluctuations. 
The NS mass and spin evolution in more complicated cases, 
like the ones observed in X-ray transients (see Campana et 
al. 1998 for a review) or the ones discussed by Burderi et al. 
(2001), are not considered here and will be discussed in a 
following paper. We considered both the case in which the 



accretion process does not excitate the r-mode instability, 
and the case in which the spin evolution is influenced by the 
gravitational waves emission. 

In figure 3 we plot the evolution of these systems in the 
gravitational mass - radius plane. The area of this plane 
that a stable NS equilibrium configuration can span is lim- 
ited from below by the sequence of stable non-rotating equi- 
librium configurations, which are the ones with minimum 
radius for a given gravitational mass, and from above by the 
mass-shedding sequence - that is the sequence of equilibrium 
configurations for which the NS angular velocity equals the 
Keplerian angular velocity at the NS surface radius. These 
two stability lines are connected by a third curve, the secu- 
lar stability curve, which limits from the right the region in 
which equilibrium configurations are stable against gravita- 
tional collapse. 

The system evolution without gravitational waves emis- 
sion are indicated by the thick solid and dashed lines in 
figure 3, where the numbered dots give the corresponding 
age and spin period. We see that the first phase of accre- 
tion (thick solid lines) is similar in all cases: a rapid spin up 
brings the star to the mass-shedding limit. Thus the star re- 
mains on the mass shedding sequence in the M-R diagram; 
it will detach from it only when the accretion stops and the 
NS switches on as a radio pulsar (thick dashed lines). 

When we consider the spin-down due to r-modes excita- 
tion, all evolutionary tracks (dot-dot-dashed lines in figure 
3) are even more similar: after a brief spin up phase, the 
NS attains the critical spin period of 2 ms. After that, the 
spin-up spin- down cycle begins and it ends either in a direct 
collapse to a black hole or in a radio pulsar phase. In both 
hypotheses, the amount of mass accreted before the pulsar 
lights up will determine the fate of the system. 

In figure 4 we show the possible behaviors of NS angu- 
lar velocity as a function of gravitational mass for given 
values of baryonic mass. From equation J]J we see that 
during the pulsar phase the spin evolution of the star will 
only depend on the baryonic mass of the star and on the 
sign of the derivative (9w/<9Mg)m b . If there is a non- 
rotating stable configuration for the given baryonic mass, 
then (du/dMG)M B > always and the pulsar will spin down 
until it stops. If such a configuration does not exists for the 
given baryonic mass, two different evolutionary tracks are 
possible: either the pulsar spins down until it comes close to 
the stability limit to gravitational collapse, where it spins up 
rapidly (in fact (du)/dMc)M B — * —oo at the instability), or 
it spins up until it collapses, with no spin down phase (i.e. 
{duo I '8Mg)m b < always). We say that in the former case 
the pulsar lights up in the supramassive regime, while in the 
latter it lights up in the extremely supramassive regime. The 
spin-up phase, that occurs in correspondence to a loss of ro- 
tational energy, is caused by the fact that, near the onset 
of the instability and in the extremely supramassive regime, 
a little loss of energy (gravitational mass) corresponds to a 
rapid decrease of the NS radius. Since the momentum of in- 
ertia of the NS strongly depends on NS radius R and on its 
compactness GM G /Rc 2 (I ~ 0.21M G R 2 /{1 - 2GM G /Rc 2 ), 
see Ravenhall & Pethick 1994), in these regimes the loss of 
rotational energy is achieved by means of a reduction of the 
momentum of inertia rather than via a spin down of the 
NS, that, instead, will spin up to partially compensate the 
reduction of the momentum of inertia. 
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Figure 3. Evolution of the NS in the binary systems of figure 1 in the gravitational mass-radius plane. The mass is in solar masses 
and the radius is in kilometers. The numbered dots indicate evolutionary stages with system age since the start of the accretion in units 
of 10 s yr (external numbers) and the spin period in milliseconds (internal, bold numbers). These numbers refer to simulations in the 
absence of gravitational waves emission. The thin dot dashed line is the maximum rotation limit, the dashed thin line indicates the 
stable nonrotating configuration and the thin solid line indicates the stability limit to gravitational collapse. The thick solid line marks 
the evolution during the accretion phase in absence of gravitational waves emission, the thick dot dot dashed line marks the evolution 
during the accretion phase, with gravitational waves driven by the r-modes instability cycle, while the thick dashed line marks the pulsar 
phase. In case a) the pulsar evolution is shown until the star reaches a spin period comparable to that of the fastest pulsar observed, 
PSR B1937+21, that is 1.56 ms and then slowly brakes to 16 ms in ~ 8 X 10 11 yrs. In case b) pulsar evolution, after a relatively short 
radio-pulsar phase lasting 5 X 10 7 yrs, brings the star to gravitational collapse as it crosses the stability limit. In case c) there is no pulsar 
phase since the accretion of matter brings the star directly to gravitational collapse. Case d) features the same companion star as case 
c), but the NS is governed by the ultrastiff EOS N. The evolution is similar to case a). Little jumps along the evolution curve are due 
to the resolution of the grid of relativistic values we used. 



In case a, in the absence of r-mode excitation, the pulsar 
lights up, with a period of 0.71 ms, when the star has a 
gravitational mass larger than the maximum non-rotating 
mass i.e. A/ Sta t = 1.8O3M0 (see figure 4i). As the pulsar loses 
energy due to dipole radiation according to equation @, it 
leaves the mass-shedding sequence and returns in the normal 
gravitational mass range, below M sta t- Thus the pulsar will 
slow down (see figure 3a), without collapsing, until it almost 
stops on a very long time-scale (it reaches 16 ms in ~ 8 x 10 1 
yr). If the r-modes are excited during the accretion process, 
the star lights up as a pulsar with the same baryonic mass 
(since the same amount of mass has been accreted), but with 
a different gravitational mass, 1.792 Mq instead of 1.805 
Mq: the larger rotational period of the newborn pulsar, 3.31 
ms, with respect to the period of 0.71 ms of the evolution 
in absence of gravitational waves, means that there is not 
enough energy to increase the gravitational mass above the 



maximum non-rotating mass. In this case, the pulsar slows 
down on a similar long timescale, reaching 16 ms in ~ 7 x 
10 11 yr. 

In case 2, if r-modes are not excited the star becomes 
a pulsar when it is in the extreme supramassive regime (see 
figure 3b) and thus equations £J imply that the NS is spin- 
ning up rather than spinning down. As it loses energy, how- 
ever, the star begins to shrink and heads towards the secular 
instability limit which brings the NS to the gravitational col- 
lapse. Its radio pulsar phase, characterized by periods well 
below 1 ms (its initial spin period is 0.581 ms, its period 
when it crosses the stability line is 0.542 ms) and by the un- 
usual sign of the period derivative, lasts ~ 5 x 10 7 yr. Such a 
short lifetime, if compared with typical lifetimes for spinning 
down pulsars, is due to the positive feedback we get in this 
case for the second of equations J7J : gravitational mass loss 
causes a spin up, which in turn will make the term on the 
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Figure 4. Different evolutionary tracks for a pulsar in the Mq — (*J plane depending on its baryonic mass, i) NS with a baryonic 
mass equal to that at the end of accretion in the evolutionary case a. The baryonic mass is low enough so that the rotating star has a 
stable nonrotating counterpart and the derivative (dui/dMo)M B is always positive, ii) NS in the supramassivc sequence. In this case, 
(Olo/OMq ) m b —* — oo near the onset of the secular instability, as noticed by various authors, but is still positive away from it. The spin 
up phase for a pulsar of this mass, assuming a magnetic field of ~ 10 8 G, lasts ~ 10 6 yr. iii) NS in the extremely supramassivc regime, 
with a baryonic mass equal to that at the end of accretion in evolutionary case b. In this case, the mass is so high that (dui / 8Mq)m b < 
for any stable configuration. The spin up phase, as can be seen from our simulations, lasts ~5x 10 7 yr. 



right of this equation bigger, causing a further mass loss. If 
the r-modes are excited no such behavior arises, since spin of 
the NS is kept above 2 ms, the maximum sustainable mass 
is smaller than that for a NS rotating at mass shedding (the 
maximum sutanaible gravitational mass for a NS spinning 
at a period of 2 ms is 1.817M©). Thus the star collapses to 
a black hole ~ 5 x 10 8 yr before the end of the accretion 
process. 

In case c, the star never lights up as a pulsar, as too 
much matter is accreted on the star and the maximum mass 
limit is exceeded. Thus the NS will directly collapse to a 
black hole (see figure 3c) . 

All this behaviors are strongly dependent on the EOS 
adopted to describe the ultradense matter. As a compari- 
son we studied the evolution of a system consisting of the 
same companion star as in case c and a \AMq NS with an 
ultra-stiff EOS (EOS N by Walecka & Serot, see Arnett and 
Bowers 1977), which has a maximum non-rotating mass of 
2.634Mq. The result of this evolution different, as we show 
in figure 3d; in this case (hereafter case d) no spinning up 
pulsar shows up, nor any accretion induced collapse hap- 
pens. Instead, we end up with a spinning down submillisec- 



ond pulsar with a very long lifetime, comparable with the 
one we obtained in the first case. If the r-modes are excited 
(see again figure 3d) the minimum attanaible period is lim- 
ited to 2 ms, and the pulsar lights up with a period of 2.27 
ms. In a following paper, we will discuss in full detail the 
effects of the EOS of the ultradense matter on the evolution 
of LMXBs. 



4 CONCLUSIONS 

We have shown that, depending on the characteristics of the 
system, especially on the amount of mass accreted on to the 
NS, on the EOS adopted to describe the Ns matter, and 
on the excitability of the r-modes, LMXBs can have quite 
different fates: they can light up as a spinning down radio 
pulsars, they can directly collapse to a black hole during the 
accretion phase or if, at the end of the accretion phase, the 
NS is left in the extreme supramassive regime, it will light up 
as an exotic, spinning up submillisecond radio pulsar with a 
relatively short lifetime. 

It is then evident that, in the hypothesis of a con- 
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servative mass transfer from the companion onto a low- 
magnetized NS and in absence of r-modes excitation, the 
accretion process, if the amount of mass accreted is not 
enough to collapse into a black hole, will end with a very 
fast spinning object, as it has been suggested before (Cook, 
Shapiro, & Teukolsky 1994a). If the NS becomes then de- 
tectable as a radio pulsar, it will have a spin period well 
below one millisecond. In fact our simulations show (cases a 
and d) that we obtain submillisecond pulsars with long life- 
times (in the former case the pulsar lifetime before reaching 
a period as long as that one of the fastest millisecond pulsar 
known to date, PSR B1937+21, is ~ 3 x 10 9 yr, while in the 
latter it becomes ~ 5 x 10 9 yr). 

If r-modes are excited by accretion, pulsars are con- 
strained to spin slower than a critical frequency, and this 
could explain why no NS spinning at submillisecond periods 
has been observed to date. However, in this situation any 
binary system in which enough mass is transferred from the 
companion to the NS will collapse to a black hole, without 
lighting up as a pulsar. Thus pulsar formation could be much 
less favoured than in other cases. 

It is likely that millisecond pulsar systems like the ones 
observed to date (i.e. systems with P > 1 ms) originate from 
different binary evolution scenarios, in which some critical 
mechanism has prevented the accretion process to continue 
until a mass as large as a significant fraction of a solar mass 
has been transferred. It is probable that the magnetic field 
of the NS has values much higher than the value we chose 
(at least at the beginning of the accretion), so that the inner 
edge of the disc could be outside the corotation radius for at 
least part of the evolution and magnetic torques could play 
an important role in the spin evolution of the NS. Moreover, 
systems in which the mass transfer rate has large fluctua- 
tions will light up as pulsars before the end of the accretion 
process, losing a large amount of mass in a so-called radio 
ejection phase as proposed by Burderi et al (2001). We will 
investigate the evolution of such systems in a future paper. 

Although there are selection effects that could have pre- 
vented the discovery of submillisecond pulsars (Burderi et al 
2001), if a self-limiting mechanism like the r-mode instabil- 
ity does not operate, submillisecond radio pulsars should 
exist and should be detectable in the future. On the other 
hand we have shown (case 3) that if the NS matter is gov- 
erned by a moderately stiff EOS like FPS (i.e. with a max- 
imum non-rotating mass of < I.9M0), the mass transfer 
can end in an accretion induced collapse to a black hole if 
as much as 1M is accreted. Although not much in known 
on the range of progenitor masses fot today's population of 
LMXB in the Galaxy, a recent study by Pfahl, Rappaport 
and Podsialowski (2003) argue that a great fraction of ob- 
served LMXBs may have descended from intermediate mass 
X-ray binaries, that is form systems with initial donor mass 
>1.5Mq. If this is true, and if we assume that in such a 
system we have equal probability that the companion trans- 
fers any amount of mass between 0.5 Mq and 2 Mq on to 
a 1.4 M Q NS, for NS governed by EOS N we have a 100% 
probability of obtaining a submillisecond pulsar, while for 
NS governed by EOS FPS we have only a 20% probabil- 
ity of obtaining a spinning down millisecond pulsar, a 9% 
probability of the formation of a spinning up submillisecond 
pulsar (doomed to gravitational collapse) and a 71% proba- 
bility of direct collapse into a black hole during the accretion 



phase. Therefore if the EOS governing NS matter is soft, a 
conservative mass transfer is more likely to end with a di- 
rect accretion-induced collapse to a black hole than with the 
formation of a submillisecond radio pulsar, and thus submil- 
lisecond pulsars could be hard to detect because of their low 
formation probability. 

On the other hand, if the r-modes are excited, the spin 
period will remain well above one millisecond for all of the 
evolution. However, making the same assumptions as before, 
the probability of forming a pulsar drops to 10% for EOS 
FPS, while for EOS N we still have a very high probability, 
~ 90%. Therefore this scenario, in which the mass transfer 
is conservative but the spin frequency is limited by the emis- 
sion of gravitational waves, implies that the EOS is stiff in 
order to have an high probability of formation of millisecond 
pulsars. 

Teherefore we should predict that the EOS of NSs is 
very stiff in order to explain the observational evidence 
(MSP are formed), if gravitational waves are indeed emit- 
ted due to r-modes excitation, while we should predict that 
the EOS of NSs is soft if they are not emitted, so that sub- 
millisecond pulsars are very uncommon, as the observations 
seem to indicate. If future observations will allow to con- 
strain the stiffness of the NS EOS on an observational basis, 
this will give an indication on wherther r-modes are indeed 
excited in LMXBs or not. 

We have shown that, if r-modes are not excited in 
LMXBs, the accretion process can leave us with an ex- 
tremely supramassive NS, that will spin up during all of 
its life as a radio pulsar (case b). It is evident that, being 
the critical baryonic mass for getting to the extremely supra- 
massive regime M cr i t an EOS-dependent feature, in principle 
the observation an accelerating (or braking) submillisecond 
pulsar can allow to exclude several EOS on an observational 
basis. 

To clarify how such an effect can help to constrain the 
EOS of ultradense matter we need to introduce the new 
concept of a critical spin period P C rit that, together with 
the minimum period (P m i n = 27r/cj max ), is peculiar to each 
EOS. Pcrit is the period below which the EOS allows only 
extremely supramassive stable configurations. In figure 5 we 
show sequences of equilibrium configurations with constant 
spin period, together with the critical baryonic mass se- 
quence. It is evident from the figure that P cr i t is equal to 
the minimum allowed period to avoid gravitational collapse 
if the star has Mb = M cr it. In fact, any constant period 
sequence P < P cr i t will only include stars of baryonic mass 
greater than the critical one. Thus any NS with P < P cr i t 
will accelerate as a consequence of energy loss due to mag- 
netic dipole radiation. Being P cr i t EOS-dependent, the de- 
tection of a submillisecond radio pulsar and the determina- 
tion of the sign of its period derivative will allow to effec- 
tively constrain the equation of state governing ultradense 
matter. 

Thus the detection of a submillisecond radio pulsar can 
impose two constraints on the EOS of the NS: 

(i) the spin period must be larger than the minimum al- 
lowed period, i.e. the spin period of the maximum rotation 
configuration, P m i n . 

(ii) if the period is shorter than P cr it, the radio pulsar 
must spin up rather than spin down. 
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Figure 5. Sequences of NSs with the same spin period are plot- 
ted as dashed lines in the gravitational mass-radius plane, while 
sequences of NSs with the same baryonic massarc plotted as solid 
lines. We consider NSs governed by EOS FPS. The dashed-dotted 
line limits stable configurations. We plot with a thick solid line 
the sequence of Ns with baryonic mass equal to the critical mass , 
M CI - lt = 2.33-Mq. Any star with Mg > M cr ; t is extremely supra- 
massive, i.e. it spins up under magnctodipolc radiation. It is inter- 
esting to note that while sequences of constant baryonic mass al- 
ways have the same shape in the gravitational mass-radius plane, 
bending from left to right with increasing radii, sequences of con- 
stant spin frequency have completely different topologies below 
the critical mass and above it (see for example the sequence with 
P = 0.7 ms and the one with P = 0.6 ms). During the radio pul- 
sar phase the star moves along a sequence of constant baryonic 
mass, decreasing its gravitational mass. It moves therefore from 
top right of the figure to the bottom left. This implies that the 
pulsar spins down as long as constant spin frequency sequences 
bend from top left to bottom right (as the one with P = 0.7 ms), 
and that it spins up if the constant spin frequency sequences it 
crosses bend from bottom left to top right in the plane (as the 
one with P = 0.6 ms does). As shown in the figure, any stable 
NS attaining a period P < P CIlt = 0.6 ms (i.e. any star who lies 
on the right of the thick dashed line) has Mg > M cr ; t . Thus any 
NS governed by EOS FPS attaining a period < 0.6 ms will spin 
up once it becomes a pulsar. 



Both P m in and P cr it are EOS dependent and are longer for 
stiffer EOSs. In fact, the detection of a submillisecond radio 
pulsar with spin period P b s undergoing a spin up will rule 
out all the stiff EOSs with Pmin > Pobs- On the other hand 
the detection of a spinning down submillisecond radio pul- 
sar, with spin period P bs, will allow us to rule out all the 
stiff EOSs with Pcrit > Pobs, because they cannot explain a 
spinning down radio pulsar with such a short spin period. 
In this case the limit is more stringent because P cr i t > P m in! 
As an example, suppose that a spinning down radio pulsar 
with a period of 0.713 ms (like the one we obtain in case 
a) will be detected: this will allow us to rule out EOS N, 
since although the minimum period for this EOS is 0.69 ms, 
any radio pulsar governed by EOS N with such a low spin 
period will spin up, being for EOS N M cr i t = 3.63M and 
Pcrit = 0.74 ms. 

In summary, in this paper we presented the first results 
obtained with a new code that allows to study in details the 
binary system evolution and the spin evolution of the NS, on 
the basis of fully relativistic calculations. We used this code 
to study the evolution of systems with conservative mass 
transfers and confirmed that the large amount of matter that 
is transferred on to the NS will spin it up to periods well 



below one millisecond, unless the emission of gravitational 
waves dissipates the excess of angular momentum. However 
in this last case the amount of mass acreted onto the NS is 
easily big enough to cause a direct collapse to a black hole. 
Therefore we concluded that presumably the recycled sys- 
tems which give origin to the MSP observed to date should 
have origin from systems with a highly nonconservative mass 
transfer. We showed that if the EOS of ultradense matter is 
not very stiff the direct collapse to a black hole is very likely 
to happen even if the r-modes are not excited. This could 
explain the lack of any observation of submillisecond radio 
pulsars even without invoking gravitational waves emission. 
As a last remark, since we showed that there is the possibil- 
ity of obtaining from binary evolution some unusual, accel- 
erating submillisecond radio pulsars, we introduced the new 
concept of the critical spin period P cr it, peculiar to each 
EOS, that can allow to effectively constrain the EOS of NS 
matter if a radio pulsar with a period below one millisecond 
will be observed in the future. 
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